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"Light bullets" are multi-dimensional solitons which are 
localized in both space and time. We show that such solitons 
exist in two- and three-dimensional self-induced-transparency 
media and that they are fully stable. Our approximate an- 
alytical calculation, backed and verified by direct numerical 
simulations, yields the multi-dimensional generalization of the 
one- dimensional Sine-Gordon soliton. 



PACS numbers: 42.50 Rh, 42.50 Si, 3.40 Kf physics/0002006 



The concept of multi-dimensional solitons that are 
localized in both space and time, alias "light bullets" 
(LBs), was pioneered by Silberberg [jlj, and has since 
then been investigated in various nonlinear optical media, 
with particular emphasis on the question of whether these 
solitons are stable or not. For a second-harmonic gener- 
ating medium, the existence of stable two- and three- 
dimensional (2D and 3D) solitons was predicted as early 
as in 1981 [gj, followed by studies of their propagation 
and stability against collapse J^-|(| , and of analogous 3D 
quantum solitons 0. In a nonlinear Schrodinger model 
both stable and unstable LBs were found S and it was 
suggested that various models describing fluid flows yield 
stable 2D spatio-temporal solitons |l . Recently, the first 
experimental observation of a quasi-2D "bullet" in a 3D 
sample was reported in Ref. ]l0| . 

In this letter we predict a new, hitherto unexplored 
type of LBs, obtainable by 2D or 3D self-induced trans- 
parency (SIT). SIT involves the solitary propagation of 
an electromagnetic pulse in a near-resonant medium, irre- 
spective of the carrier-frequency detuning from resonance 
fjl|Q. Tne SIT soliton in ID near-resonant media (TJ] is 
exponentially localized and stable. In order to investigate 
the existence of "light bullets" in SIT, i.e. solitons that 
are localized in both space and time, one has to consider 
a 2D or 3D near-resonant medium. Here we present an 
approximate analytical solution of this problem, which 
is checked by and in very good agreement with direct 
numerical simulations. 

Our starting point are the two-dimensional SIT equa- 
tions in dimensionless form Il4| 



l£xx ~!~ £% F 

V T - £ W = 
W T + ~{£*V + V*£) = 0. 



(la) 
(lb) 

(lc) 



Here £ and V denote the slowly-varying amplitudes 
of the electric field and polarization, respectively, W is 
the inversion, z and x are respectively the longitudinal 
and transverse coordinates (in units of the effective ab- 
sorption length a c ff), and r the retarded time (in units 
of the input pulse duration t p ). The Fresnel number F 
(F > 0), which governs the transverse diffraction in 2D 
and 3D propagation, is incorporated in x and the detun- 
ing AJ7 of the carrier frequency from the central atomic 
resonance frequency is absorbed in £ and V [^5|. We 
have neglected polarization dephasing and inversion de- 
cay, considering pulse durations that are much shorter 
than the corresponding relaxation times. Eqs. (Q) are 
then compatible with the local constraint \P\ 2 + W 2 = 1, 
which corresponds to conservation of the Bloch vector 

& 

The first nontrivial question is to find a Lagrangian 
representation for these 2D equations, which is necessary 
for adequate understanding of the dynamics. To this end, 
we rewrite the equations in a different form, introducing 
the complex variable <fi defined as follows J16|| 



l + W 



p* 



1 - W 



V 



20* 



W 



h 1 
(2) 



Eqs. (lb) and fllq ) can then be expressed as a single equa- 
tion, cj) T + {£/2)4> 2 + (l/2)£* = 0. Next, we define a vari- 
able / so that 4> = 2/ T / (£f). In terms of /, the previous 
equation becomes f TT - (£ T /£)f T + (1/4) \£\ 2 f = 0. This 
equation is equivalent to 



fr = - 2 £g 



(3a) 
(3b) 



with g = / (f>. Applying the same transformations to 
Eq. (|lj) yields 



i£xx + £z- 2fg* = 0. 



(4) 



The Lagrangian density corresponding to Eqs. ([|) and 
(Q) can now be found in an explicit form, 

C(x, t) - - A £ x £* x + '-{££; - £ z £*) - % - (f*g£ - fg*£*) 



ff* - ff 



{99 -99 ) 



(5) 
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Now we proceed to search for LB solutions. Before 
resorting to direct simulations, we obtain an analytical 
approximation of the solutions. The starting point for 
this approximation is the well-known soliton solution for 
ID SIT (the Sine-Gordon soliton) 



£(t, z) — ±2 asechO 
V(t,z) = ±2sech6tanh6 
W(t, z) = sech 2 6 - tanh 2 9, 



(6a) 
(6b) 
(6c) 



with 9(r, z) — ar — ^ + 9o, and a, 9o arbitrary real 
parameters. Equation (6a) is also called a 27r-pulse, be- 
cause its area £(t, z)dr = ±2ir. 

Returning to the 2D SIT equations, we notice by 
straightforward substitution into Eqs. (||) that a 2D so- 
lution with separated variables, in the form £(t,z,x) = 
£i{t, z) £2(2;) (and similarly for / and 5), does not exist. 
To look for less obvious solutions, we first split equa- 
tions (|l|) into their real and imaginary parts, writing 
£ = £ x + i£ 2 and V = Vi + iV 2 - 



£ixx + £\z — V\ — 
£ixx — £iz + P2 = 

V Xt - £ 1 W = 

V 2 r ~ £2 W = 
W T + £{Px + £ 2 V 2 = 0. 



(7a) 
(7b) 
(7c) 
(7d) 
(7e) 



In the absence of the x-dependence, these equations are 
invariant under the transformation (£\, Vx) <-* (£2,^2) ■ 
This suggests a ID solution in which real and imaginary 
parts of the field and polarization are equal, £ 1 = £2 and 
V\ = V21 and such that the total field and polarization 
reduce to the SG solution (||). Our central result is an 
approximate but quite accurate (see below) extension of 
this solution, applicable to the 2D SIT equations. In 
terms of the original physical variables it is given by 



£{t,z,x) 
T>(t,z,x) 



±2a \J sechOisech02 exp(— iAQr + iir/A) 
± ^/seche 1 sech6 2 { (tanh 9 1 + tanh9 2 ) 2 + 

ia 2 C 4 [(tanh9i - tanh9 2 ) 2 - 



2(sech 2 9i + sech 2 9 2 )] 2 } 1/2 exp(-iA^r 



W(t,z,x) = [1 — sech9isech92{(tanh 9i - 
ia 2 C 4 [(tanh9i - tanh9 2 ) 2 



tanh9 2 ) 5 



2(sech 2 9i +sech 2 9 2 )] 2 }] 



1/2 



with 



z 

91 = ar 

a 

92 = ar 

a 



9 + 

9 - 
arctan (J 3 2/ 'Pi) 



Cx 
Cx, 



sat- 



Here a, 9o and C are real constants. Equations y£ 
isfy the two-dimensional SIT equations (|7a|) and ( |7b[ ) and 



obey the normalization condition V\ + V\ + W 2 = 1. 
They reduce to the Sine-Gordon solution for C = 0. 
The accuracy to which Eqs. (j|) satisfy Eqs. (|7^)-(|7e]) is 
0(aC 2 ), which requires that \a\C 2 < 1. This is the 
single approximation made. Numerical simulations dis- 
cussed later on verify that Eq. (j|) indeed approximates 
the exact solution of Eq. (0) to a high accuracy. In ad- 
dition, we have checked that substitution of (||) into the 
Lagrangian ^ and varying the resulting expression with 
respect to the parameters a and C yields zero. This 
"variational approach" is commonly used to obtain an 
approximate "ansatz" solution to a set of partial differ- 
ential equations in Lagrangian representation [jl8| . Equa- 
tions (||) represents a light bullet, which decays both in 
space and time and is stable for all values of z. The latter 
follows directly from ( |Sa|) and also from the Vakhitov- 
Kolokolov stability criterion |19| . 




FIG. 1. The electric field in the 2D "light bullet", \£\, as a 
function of time r (in units of the input pulse duration t p ) and 
transverse coordinate x (in units of the effective absorption 
length a e ff) after propagating the distance z — 1000. Param- 
eters used correspond to a = 1, C = 0.1 and Go = 1000. 



(8a) 



(8c) 
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FIG. 2. Contourplot of Fig. |l] in the (r,a;)-plane. Regions 
with lighter shading correspond to higher values of the electric 
field. Note the different time scale than that of Fig. pi. 




FIG. 3. The polarization in the 2D "bullet", \V\, as a func- 
tion of time r and transverse coordinate x. Parameters used 
are the same as in Fig. |l| 

Figs. show the electric field and polarization, gener- 
ated by direct numerical simulation of the 2D SIT equa- 
tions ([!]) at the point z = 1000, using (||) as an initial 
ansatz for z = 0. To a very good accuracy (with a de- 
viation < 1%), they still coincide with the initial con- 
figuration and analytic prediction (||). The electric field 
has a typical shape of a 2D LB, localized in time and 
the transverse coordinate x, with an amplitude 2a and 
a nearly sech-form cross-section in a plane in which two 
of the three coordinates r, z and x are constant. The 
ratio C/a determines how fast the field decays in the 
transverse direction. For \C/a\ <C 1 (then \C\ < 1, as 
\a\C 2 -C 1), we have a relatively rapid decay in r and 
slow fall-off in the ^-direction, as is seen in Fig. [|. In the 
opposite case, \C/a\ 3> 1, the field decays more slowly in 
time and faster in x. The polarization field has the shape 
of a double-peaked bullet. Its cross-section at constant x 
displays a minimum at mm ~ 0, where |^(0 m in)l ~ 0, 
and maxima at 9± = ±Arcosh(v / 2), where \V(Q±)\ « 1. 
The field and polarization decay in a similar way, which 
is a characteristic property of SIT |Q. Also the inver- 
sion decays both in time and in x, but to a value of — 1 
instead of zero, corresponding to the atoms in the ground 
state at infinity. A numerical calculation of the field area 
at x = yields dr\£ (r, z, 0)| = 6.28 ± 0.05 « 2tt, 
irrespective of z. By analogy with the SG soliton, one 
might thus name this a " 2-7T bullet" . 

We have also numerically obtained axisymmetric stable 
LBs in a 3D SIT medium, see Fig. ||. The 3D medium is 
described by Eqs. (P with the first one replaced by 

-i(£ rr + r- l E r )+E z -V = 0, (9) 

where r = y/x 2 + y 2 is the transverse radial coordinate. 
Searching for an analytic 3D bullet solution in the trans- 



verse plane proves to be difficult. However, in the limit 
of either large or small r, an approximate analytic so- 
lution may be found. For large r, it again takes the 
form (0), but now with ©i = ar — z/a + 0rj + Cr and 
02 = ar — z/a + ©o — Cr, where a, 0o, and C are 
constants, \a\C 2 <C 1, and it is implied r 3> 1/|C|. It 
is in sufficiently good agreement (deviations < 5%) with 
results of simulation of the 3D equations, using this so- 
lution as an initial ansatz. Comparison of Figs. [I] and |] 
shows that the 2D and 3D bullets have similar shapes, 
but the 3D one decays faster in the radial direction for 
small r than the 2D bullet in its transverse direction. 



E 




FIG. 4. The electric field in the 3D "light bullet", \£\, as a 
function of time r and transverse radial variable r after prop- 
agating the distance z = 1000. Parameters used correspond 
to a = 1, C = 0.1 and 9 = 1000. 

For constant r, the 2D and 3D bullets are localized 
in both the propagation direction z and the transverse 
direction(s). One may also ask whether there exist SIT 
solitons which are traveling (plane) waves in z and lo- 
calized in x (and y). Using a symmetry argument, it is 
straightforward to prove that they do not exist. Starting 
from the SIT equations ([!]) (in 2D, the 3D case can be 
considered analogously) we adopt a plane- wave ansatz for 
E and V, changing variables as follows: x — ► \fkx (assum- 
ing k > 0), £(t,z,x) — > £(t,x) exp(— ikz), V(t,z,x) — > 
k^ 1 V(t 7 x) exp(— ikz), and W(t,z,x) — * k~ x W(t,x). 
The equations for the real and imaginary parts of the 
field then become 

£2 XX -£ 2 -Vi = (10a) 

£xxx - £i + V 2 = o, (10b) 

with the equations for V T and W T given by (^q)-(|7e|). Us- 
ing the transformation (£ 1 ,'Pi) «-» (£2^2)1 which leaves 
the last three equations invariant but changes the first 
two, one immediately finds that (^|) only admits the 
trivial solution £ x = £ 2 = Vi = "P 2 = 0, W = -1. 

The observation of " light bullets" in a SIT process re- 
quires high input power of the incident pulse and high 
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density of the two-level atoms in the medium, in order 
to achieve pulse durations short compared to decoher- 
ence and loss times. These requirements are met e.g. 
for alkali gas media, with typical atomic densities of 
~ 10 11 atoms/cm 3 and relaxation times ~ 50 ns [ p0[ , 
and for optical pulses generated by a laser with pulse 
duration r p < 0.1 ns. In order to include transverse 
diffraction, the incident pulse should be of uniform trans- 
verse intensity and satisfy a e ed 2 /X < 1 p0| , where A and 
d are its carrier wavelength and diameter respectively 
p0{ . The parameter a in the solution (||), which deter- 
mines the amplitude of the bullet and its decay in time, 
corresponds to a ~ k z t p v p [T9], with k z the wavevec- 
tor component along the propagation direction z and v p 
the velocity of the pulse in the medium, and can thus 
be controlled by the incident pulse duration and veloc- 
ity. The parameter C ~ n x L x , where n x is the trans- 
verse component of the wavevector and L x is the spa- 
tial transverse width of the pulse, is also controlled by 
the characteristics of the incident pulse and should sat- 
isfy the condition k z k x L z L x <C 1. For a homogeneous 
(atomic beam) absorber, the effective absorption length 
a c s ~ 10 4 m _1 and the Fresnel number F can range from 
1 to 100 p0]| . The bullets then decay on a time scale of 
t ~ 1 — 10 t p ~ 10 ns and transverse length of x ~ 0.1 — 1 
mm, which is well within experimental reach. 

In conclusion, we predict the existence of fully stable 
"light bullets" in 2D and 3D self- induced transparency 
media. The prediction is based on an approximate 
analytical solution of the multi-dimensional SIT equa- 
tions and verified by direct numerical simulation of these 
PDE's. Our results suggest an experiment aimed at de- 
tection of this "bullet" in an SIT- medium and opens the 
road for analogous searches for "light bullets" in other 
nonlinear optical processes, such as, e.g., stimulated Ra- 
man scattering, which is analogous to SIT. 
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